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Abstract 



A fractal construction shows that, for any (3 > 0, the /3-skeleton of a point set can have arbitrarily 
large dilation. In particular this applies to the Gabriel graph. 

1 Introduction 

A number of authors have studied questions of the dilation of various geometric graphs, defined as the 
maximum ratio between shortest path length and Euclidean distance. 



o 

CN 

o 

u 

For instance, Chew [||] showed that the rectilinear Delaunay triangulation has dilation at most VTO and 
that by placing points around the unit circle, one could find examples for which the Euclidean Delaunay 
triangulation has dilation arbitrarily close to w/2. In the journal version of his paper [Q], Chew added 
a further result, that the graph obtained by Delaunay triangulation for a convex distance function based 
on an equilateral triangle has dilation at most 2. Chew's conjecture that the Euclidean Delaunay dilation 
was constant was proved by Dobkin et al. [§], who showed that the Delaunay triangulation has dilation at 
£ — . 1 most ipir where ip is the golden ratio (1 + y/5)/2. Keil and Gutwin [ |10| ] further improved this bound to 

Q ■ , 2 7 « 2.42. 

Q\ . 3cos(tt/6) 

, Das and Joseph [|4j] showed that these constant dilation bounds hold for a wide variety of planar graph 

""^ | construction algorithms, satisfying the following two simple conditions: 

O ' 

• Diamond property. There is some angle a < it, such that for any edge e in a graph constructed 
by the algorithm, one of the two isosceles triangles with e as a base and with apex angle a contains 
no other site. This property gets its name because the two triangles together form a diamond shape, 
depicted in Figure [j](a). 

• Good polygon property. There is some constant d such that for each face / of a graph constructed by 
the algorithm, and any two sites u, v that are visible to each other across the face, one of the two paths 
around / from u to v has dilation at most d. Figure [j](b) depicts a graph violating the good polygon 
property. 

Intuitively, if one tries to connect two vertices by a path in a graph that passes near the straight line 
segment between the two, there are two natural types of obstacle one encounters. The line segment one is 
following may cross an edge of the graph, or a face of the graph; in either case the path must go around 
these obstacles. The two properties above imply that neither type of detour can force the dilation of the pair 
of vertices to be high. 

*Dept. of Information and C omputer Science, University of California, Irvine, CA 92697-3425, eppstein@ics.uci.edu, 
http://www.ics.uci.edu/^eppstein^ 
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Figure 1 . (a) Diamond property: one of two isosceles triangles on edge is empty, (b) Graph violating good polygon 
property: ratio of diagonal to boundary path is high. 



For a survey of further results on dilation, see Our interest here is in another geometric graph, the (3- 
skeletons [11, 13], which have been of recent interest for their use in finding edges guaranteed to take part in 
the minimum weight triangulation [p], Q [l4|] As a special case, (3 = 1 gives the Gabriel graph, a subgraph of 
the Delaunay triangulation and the relative neighborhood graph, and a supergraph of the minimum spanning 
tree. These graphs have a definition (given below) closely related to Das and Joseph's diamond property. The 
value (3 is a parameter that can be taken arbitrarily close to zero; for any point set, as beta approaches zero, 
more and more edges are added to the /3-skeleton until eventually one forms the complete graph. Therefore 
it seems reasonable to guess that, for sufficiently small (3, the /3-skeleton should have bounded dilation. 
Such a result would also fit well with Kirkpatrick and Radke's motivation for introducing /^-skeletons in 
the study of "empirical networks": problems such as modeling the probability of the existence of a road 
between cities [11]. 

In this paper, we show that this is surprisingly not the case. For any (3, we find point sets for which 
the /^-skeleton has arbitrarily high dilation. Our construction uses fractal curves closely related to the Koch 
snowflake. We show that the point set can be chosen in such a way that the /^-skeleton forms a path with 
this fractal shape; the fact that the curve has a fractal dimension greater than one then implies that the graph 
shortest path between its endpoints has unbounded length. 



2 Beta-skeletons 

The /^-skeleton [ 11 , |l3|] of a set of points is a graph, defined to contain exactly those edges ab such that no 
point c forms an angle acb greater than sin -1 1//3 (if (3 > 1) or ir — sin -1 (if /3 < 1). 

Equivalently, if f3 > 1 , the /3-skeleton can be defined in terms of the union U of two circles, each having 
ab as a chord and having diameter (3 d(a, b). Edge ab is included in this graph exactly when U contains no 
points other than a and b. 

If (3 = 1, an edge ab is included in the /3-skeleton exactly when the circle having ab as diameter contains 
no points other than a and b. The 1 -skeleton is also known as the Gabriel graph [§]. 

If < (3 < 1, there is a similar definition in terms of the intersection / of two circles, each having ab 
as a chord and having diameter d(a, b)/(3. Edge ab is included in the /3-skeleton exactly when / contains no 
points other than a and b. 

Figure ^ depicts these regions for (3 = V2 (union of circles), (3 = 1 (single circle), and (3=1/ \fl 
(intersection of circles). 

As noted above, /3-skeletons were originally introduced for analyzing empirical networks. Gabriel 
graphs and /3-skeletons have many other applications in computational morphology (combinatorial meth- 
ods of representating shapes). Gabriel graphs can also be used to construct minimum spanning trees, since 



2 



Figure 2. Empty regions for %/2-skeleton, Gabriel graph, and 1 / v^-skeleton 





Figure 3. Fractal curves P(yr/4, k) for k — 1 , 2, 3. 



the gabriel graph contains the MST as a subgraph. More recently, various researchers have shown that 
/3-skeletons (for certain values of > 1) form subgraphs of the minimum weight triangulation [|l], [| pfl . 

Su and Chang [12] have described a generalization of Gabriel graphs, the ^-Gabriel graphs, in which 
an edge is present if its diameter circle contains at most k — 1 other points. One can similarly generalize 
/3-skeletons to &-/?-skeletons. Our results can be made to hold as well for these generalizations as for the 
original graph classes. 



3 Fractals and dilation 



Our construction showing that beta-skeletons have unbounded dilation consists of a fractal curve with a 
recursive definition similar to that of a Koch snowflake. For a given angle 6 define the polygonal path 
P(6, 1), by following a path of five equal-length line segments: one horizontal, one at angle 9, a second 
horizontal, a segment at angle —6, and a third horizontal. 

We then more generally define the graph P(6, k) to be a path of 5^ line segments, formed by replacing the 
five segments of P(6, 1) with congruent copies of P(9,k— 1), scaled so that the two endpoints of the path are 
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Figure 4. Fractal curve is contained in a diamond. 



at distance one from each other. Figure ^] shows three levels of this construction. In the drawing of Figure [|, 
the orientations of the five copies of P(9, k— 1) alternate along the overall path, so that the horizontal copies 
are in the same orientation as the overall path and the other two copies are close to upside-down, but this 
choice of orientation is not essential to our construction. 

Note that, if we denote the length of P(9, k) by £ k = 4(0), then l x > 1 and £ k = t\. 

Lemma 1. P(8, k) is contained within a diamond shape having the endpoints of the path as its diagonal, 
and with angle 6 at those two corners of the diamond. 

Proof: This follows by induction, as shown in Figure ||, since the five such diamonds containing the five 
copies of P(9, k — 1) fit within the larger diamond defined by the Lemma. □ 

Lemma 2. IfO<(ir — sin -1 j3) /2, P(6, k) is the (i-skeleton of its vertices. 

Proof: We show that, if a and b are non-adjacent vertices in the path, then there is some c forming an angle 
of at least tt — sin -1 j3. We can assume that a and b are in different copies of P(9, k — 1), since otherwise the 
result would hold by induction. But no matter where one places two points in different copies of the small 
diamonds containing the copies of P(6, k — 1) (depicted in Figure ||), we can choose one of the three interior 
vertices of P(9, 1) as the third point c forming an angle acb > n — 29. The result follows from the assumed 
inequality relating 9 to j3. □ 

For instance, the graphs P(ir/4,k) depicted in Figure ^] are Gabriel graphs of their vertices. A more 
careful analysis shows that larger values of 9 still result in a /3-skeleton: if the orientations of the copies of 
P(9,k— 1) that form P(9, k) are chosen carefully, P(9,k) is contained in only half the diamond of Lemma [T], 
and angle acb in the proof above can be shown to be > tt — 39/2. 

Theorem 1. For any (3 > there is a c > such that ^-skeletons ofn-point sets have dilation Q(n c ). 

Proof: We have seen that we can choose a 9 such that the graphs P(9, k) are /3-skeletons. Since the endpoints 
of the path are at distance one from each other, the dilation of P(9,k) is 4 = t\. Each such graph has 
n = 5 k + 1 vertices and dilation l\ = n lo ^ e i-°W. Since l\ > 1, log 5 iy > 0. □ 

4 Upper Bounds 

We have shown a lower bound of Q,{n c ) for the dilation of /3-skeletons, where c is a constant depending on 
f3, and approaching zero as (3 approaches zero. This behavior of having length a fractional power of n is 
characteristic of fractal curves; is it inherent in /3-skeletons or an artifact of our fractal construction? We 
now show the former by proving an upper bound on dilation of the same form. 
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Figure 5. Tree of triangles and corresponding abstract tree. 



To do this, we define an algorithm for finding short paths in /3-skeletons. As a first start towards such an 
algorithm, we use the following simple recursion: to find a path from s to t, test whether edge st exists in 
the /3-skeleton. If so, use that edge as path. If not, some r forms a large angle srt; concatenate the results of 
recursively finding paths from s to r and r to t. 

For (3 < 1, sr and rt are shorter than st, so this algorithm always terminates; we assume throughout 
the rest of the section that (3 < 1. We can represent the path it finds as a tree of triangles, all having an 
angle of at least ir — sin -1 [3, rooted at triangle srt (Figure ||). The hypotenuse of each triangle in this tree is 
equal to one of the two shorter sides of its parent. Note that the triangles may overlap geometrically, or even 
coincide; we avoid complications arising from this possibility by only using the figure's combinatorial tree 
structure. We will bound the length of the path found by this algorithm by manipulating trees of this form. 
For any similarly defined tree of triangles, we define the boundary length of the tree to be the following 
formula: 

\T\ = dist(^, t) + 2_j (perim(A) — 2 • hypotenuse(A)). 

Aer 

In other words, we sum the lengths of all the short sides of the triangles, and subtract the lengths of all 
non-root hypotenuses. If the tree forms a non-self-intersecting polygon, such as the one shown in the figure, 
this is distance from s to t "the long way" around the polygon's perimeter 

Lemma 3. For the tree defined by the algorithm above, \T\ is the length of the path constructed by the 
algorithm. 

Proof: This can be shown by induction using the fact that the path from s to t is formed by concatenating 
the paths from s to r and r to t. □ 

Our bound will depend on the number of leaves in the tree produced above. However, this number may 
be very large, larger than n, because the same vertex of our input point set may be involved in triangles in 
many unrelated parts of the tree. Our first step is to prune the tree to produce one that still corresponds in a 
sense to a path in the /^-skeleton, but with a good bound on the number of leaves. 

Lemma 4. For any (3 < 1, we can find a tree like the one described above, with at most In leaves, for 
which \T\ is the length of some path in the (3-skeleton from s to t. 

Proof: Define a "leaf vertex" to be the vertex opposite the hypotenuse of a leaf triangle in T. We prune 
the tree one step at a time until each vertex appears at most twice as a leaf vertex. At each step, the path 
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Figure 6. Making single-leaf tree longer: subdivide triangles, increase angles, add children. 



Figure 7. Logarithmic spiral formed by keeping fixed angle to destination point. 

corresponding to T (and with length at most \T\) will visit all the leaf vertices in tree order (as well as 
possibly visiting some other vertices coming from interior nodes of the tree). 

Suppose some vertex v appears three or more times. Then we prune T by removing all subtrees descend- 
ing from the path between the first and last appearance of v (occurring between the two appearences in tree 
order), and we shorten the corresponding path by removing the portion of it between these two appearances 
of v. At each step, the change to \ T\ comes from subtracting some triangle short side lengths corresponding 
to the subtrees removed from T, as well as adding some hypotenuses of triangles from the same subtrees. 
Each subtracted side length that is not cancelled by an added hypotenuse corresponds to one of the edges 
removed from the path, so the total reduction in \T\ is at most as great as the total reduction in the length 
of the path, and the invariant that \T\ bounds the path length is maintained. After this pruning, there will 
be no leaves between the two appearances of v, and no new leaves are created elsewhere in the tree, so the 
invariant that the path visits the leaf vertices in order is also maintained. 

This pruning process removes at least one appearance of v, and so can be repeated at most finitely many 
times before terminating. □ 

We use induction on the number of leaves to prove bounds on \T\. The following lemma forms the base 
case: 

Lemma 5. Let T be a tree of triangles, all having an angle of at least 6 > ir/2 opposite the edge connecting 
to the parent in the tree, with exactly one leaf triangle, and scaled so that the hypotenuse of the root triangle 
has length 1. Then \T\ < — l/cos#. 

Proof: Since \T\ does not depend on the ordering of tree nodes, we can assume without loss of generality 
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Figure 8. Exponent in the bound of Theorem g, as a function of (3, 

that each node's child is on the left. For any such tree, we can increase \T\ by performing a sequence of 
the following steps: (1) If any triangle has an angle greater than 9, change it to one having an angle exactly 
equal to 9, without changing any other triangle shapes. (2) If any triangle has a ratio of left to right side 
lengths less than some value C, split it into two triangles by adding a vertex on the right side. (3) Add a 
child to the leaf of T. These steps are depicted in Figure |6[ 

The result of this sequence of transformations is the concatenation of many triangles with angles equal 
to 9, very short left sides, and right sides with length close to that of the hypotenuse. In the limit we get 
a curve from s to t formed by moving in a direction forming an angle it — 9 to t, namely the logarithmic 
spiral (Figure [7|). Integrating the distance traveled on this spiral against the amount by which the distance 
to t is reduced shows that it has the length formula claimed in the lemma. Since we reach this limit by a 
monotonically increasing sequence of tree lengths, starting with any finite one-leaf tree, any finite tree must 
have length less than this limit. □ 

More generally, we have the following result. 

Lemma 6. Let T be a tree of triangles, all having an angle of at least 9 > tt/2 opposite the edge connecting 
to the parent in the tree, with k leaf triangles, and scaled so that the hypotenuse of the root triangle has 
length 1. Then \T\ < (-l/cos0) 1+ L lo 62*J. 

Proof: We prove the result by induction on k; Lemma |5] forms the base case. If there is more than one leaf 
in T, form a smaller tree T' by removing from T each path from a leaf to the nearest ancestor with more than 
one child. These paths are disjoint, and each such removal replaces a subtree with one leaf by the edge at 
the root of the subtree, so using Lemma || again shows that \T\ < —\T'\/ cos 9. Each leaf in T' has two leaf 
descendants in T, so the number of leaves in T is at most k/2 and the result follows. □ 

This, finally, provides a bound on /3-skeleton dilation. 



Theorem 2. For (3 < \/3/2 0.866025, any ^-skeleton has dilation 0(n c 
depending on (3 and going to zero in the limit as f3 goes to zero. 



where c < I is a constant 



Proof: We have seen (Lemma ||) that we can connect any pair of vertices in the skeleton by a path with 
length bounded by \T\, where T is a tree of triangles in which all angles are at least it • '" ■ ~ 1 
T has at most 2n leaves. By Lemma 0, the length of such a tree is at most 



sin f3, and where 
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which has the form specified in the statement of the theorem. □ 

Figure || shows the growth of the exponent c as a function of (3. For v3/2 < /3 < 1, the theorem does 
not give the best bounds; a bound of n — 1 on dilation can be proven using the fact that the skeleton contains 
the minimum spanning tree. 
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